Modelling cytoskeletal traffic: an interplay between passive diffusion and active transport 
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We introduce the totally asymmetric exclusion process with Langmuir kinetics (TASEP-LK) on a network as 
a microscopic model for active motor protein transport on the cytoskeleton, immersed in the diffusive cytoplasm. 
We discuss how the interplay between active transport along a network and infinite diffusion in a bulk reservoir 
leads to a heterogeneous matter distribution on various scales. We find three regimes for steady state transport, 
corresponding to the scale of the network, of individual segments or local to sites. At low exchange rates strong 
density heterogeneities develop between different segments in the network. In this regime one has to consider 
the topological complexity of the whole network to describe transport. In contrast, at moderate exchange rates 
the transport through the network decouples, and the physics is determined by single segments and the local 
topology. At last, for very high exchange rates the homogeneous Langmuir process dominates the stationary 
state. We introduce effective rate diagrams for the network to identify these different regimes. Based on this 
method we develop an intuitive but generic picture of how the stationary state of excluded volume processes on 
complex networks can be understood in terms of the single- segment phase diagram. 

PACS numbers: 87.16.A-, 87.16.Uv, 05.60.Cd, 89.75.Hc 



Introduction Statistical physics has been very success- 
ful in deducing macroscopic properties of materials from the 
interactions between their microscopic components. Active 
matter systems on the other hand, such as active colloids, an- 
imal flocks or cytoskeletal assemblies, are prone to develop 
out-of equilibrium patterns. These spatial heterogeneities are 
in fact essential to life processes - 
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FIG. 1. (color online) Statistical physics model of cytoskeletal trans- 
port, capturing the competition between active transport of particles 
through a network (cytoskeleton) with diffusion in a bulk reservoir 
(cytoplasm). 



Here we would like to initiate a microscopic statistical 
physics approach to describe collective organization of molec- 
ular motors in cells. Motor proteins navigate actively through- 
out the cell MAM along the cytoskeleton, a network of filamen- 
tous assemblies spanning the cytoplasm. These proteins can 
exert forces, depolymerize filaments and transport biological 
cargos along the cytoskeleton [4], and thus play an important 
role in the assembly, self-organization and functioning of cells 
10]. Single-molecule properties of the motors are well-studied 
and can now be measured accurately L6] . But understanding 
how motors collectively self-organize remains a very impor- 



tant step in developing a microscopic vision of intracellular 
organization. 

Macroscopic approaches to study intracellular motor pro- 
tein transport have been developed |7|], including efforts to in- 
troduce microscopic aspects [8, 9], but a generic microscopic 
picture of cytoskeletal transport has yet to emerge. Here 
we present a tentative approach using a minimalistic model 
for cytoskeletal active transport, which consists of directed 
motion of particles (motor proteins) along the network (cy- 
toskeleton) and diffusion in the bulk (cytoplasm), see Fig. Q] 
We model the directed motion along the cytoskeleton using 
the totally asymmetric simple exclusion process (TASEP) it 1 OTI 
along a disordered directed network II ill . The binding and 
unbinding of particles between the network and the bulk is 
represented via Langmuir kinetics (LK) |[l2ll . We consider the 
particles in the bulk to be infinitely diffusive which, as we will 
show, is a relevant limiting case. 

TASEP is a fundamental model in non-equilibrium physics 
11311 . but is also used in more applied topics, such as mod- 
elling macromolecules moving through capillary vessels [14], 
electrons hopping through a quantum-dot chain II 1 311 and ve- 
hicular traffic [16]. LK on the other hand is a well known 
fundamental equilibrium process in chemical physics 112]. 

Our model constitutes a generalization of transport through 
closed networks 111 ill to open systems, as they are relevant for 
cytoskeletal transport. It may also be seen as generalization to 
a large scale network of the totally asymmetric exclusion pro- 
cess with Langmuir kinetics (TASEP-LK) on a single segment 
fnll . which has been shown to quantitatively describe in-vitro 
experiments of motor proteins H 1 811 . Our study differs from 
previous work [8] in that we consider exclusion interactions 
and disordered networks, both essential points to the physical 
picture we develop. 

The fundamental question we address is how spatial hetero- 
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FIG. 2. (color online) Sketch of the method to study transport 
through complex networks (a), here for TASEP-LK. Each segment 
of the network (b) is considered to connect to two reservoirs with ef- 
fective rates a e ff and /3 e ff (c), which are set by the junction densities. 



geneities emerge in such open active systems, here due to the 
competition between diffusion in a reservoir (which spreads 
particles homogeneously) and active transport along a net- 
work (which generates heterogeneities Willi ). 

We will show how this competition is governed by the hop- 
ping rate for particles on the network, the exchange rates with 
the reservoir and the filament length. One major result of 
our work is that there exist three regimes of transport through 
complex networks which are linked to the scale at which spa- 
tial density heterogeneities arise: the network, the segments or 
the sites. 

Microscopic model for cytoskeletal transport We repre- 
sent the cytoskeleton as a network of directed segments of L 
sites each, connected by junction sites, see Fig. [2] (a). We use 
random networks to reflect the topological complexity of the 
cytoskeleton, a standard approach in modelling networks Il9ll . 
Specifically we use Erdos-Renyi graphs of mean connectivity 
c, as in 111 ill . Whereas the specific topology is not important 
for our qualitative results, the fact that networks are irregular, 
i.e. that the number of incoming c l v and outgoing c° segments 
of the junction v differ, is relevant. 

In each directed segment particles move according to 
TASEP-LK rules iQ: particles hop uni-directionally at rate 
p, subject to exclusion interactions. Furthermore, particles 
obey binding/unbinding kinetics with attachment rate uja and 
detachment rate ujd at every site along the network. Particles 
in the reservoir are assumed to diffuse infinitely fast. 

The phase diagram of TASEP-LK has been determined Il7ll 
for a single segment connecting a reservoir with entry rate a to 
a reservoir with exit rate /3 (see supplemental materials [20]). 
TASEP-LK is best characterized in terms of the dimension- 
less parameters ft a = ^aL/p and VLd = ujdL/p, where fin 
relates the distance an isolated particle typically moves before 
detaching to the filament length 111 711 . It is convenient to con- 
sider the parameters Q = (Qa + ^d)/2, characterizing the 
total exchange between reservoir and segment, and the ratio 
K = Qa/Qd, which sets the equilibrium Langmuir density 
pi = K/(K+l). 

The (a, 0) -phase diagram of TASEP-LK is fully charac- 
terized in terms of the parameters Q and K (see 
and Fig. [4]). In the following we refer to low density (LD), 
high density (HD) and maximum current (MC, or M in Il7ll ) 



phases. All of these reduce to TASEP phases with a flat den- 
sity profile for Q = fl 1 011 . Furthermore LD and HD zones 
can coexist on the same segment separated by a domain wall 
leading to a coexistence phase (LD-HD). 

Mean field method for TASEP-LK on networks We ana- 
lyze TASEP-LK by extending the mean field arguments for 
TASEP presented in 111 ill . In this approach every segment 
(v,v f ) connecting two vertices v and v' is considered to be 
governed by effective entry and exit rates v ,\ and v ,y 
see Fig. [2] (c). These rates are in turn determined by the 
average densities p v and p v > at the junction sites 

a l S v> = P Pv/ C v and Pt% = P (1 - Pv)- Balancing the 
currents at the junctions leads to the following closed set of 
equations in p v : 



dpv 
dt 



Pv' 

c°, 

1 1 



1 - Pv 



Pv 



1 - p v , (1) 



where the sums are over incoming (outgoing) segments, and 
J ± [a/p,l3/p] are the currents entering (leaving) a segment 
with rates a (/?.) The expressions for J ± are readily available 
InlS]- Due to the Langmuir process the current is not con- 
stant along the segment such that J~ ^ J + . Solving Eqs. (Q]) 
yields the complete stationary state of all segments in the net- 
work. 

The overall particle density on a network immersed in a 
reservoir is equal to the Langmuir density pi, set by the ratio 
K ll23ll . In the following we discuss the physical phenomena 
at fixed K as the exchange parameter Vt smoothly interpolates 
from purely active transport in standard TASEP (Vt = 0) to 
the diffusion dominated limit (Q — >• 00). The data shown here 
has been obtained using K = 1.5, which is a reasonable value 
for motor proteins [ 22] and theoretically not specific in any 
way: our analysis is general, and only very high and very low 
values of K require an additional discussion |230. 

Decoupling due to particle exchange In principle the con- 
tinuity equations 0} couples the densities p v and p v > of those 
junctions which are linked by a segment (v, v')\ it is this cou- 
pling which makes the transport problem global, requiring to 
analyze the whole network simultaneously. 

Here we exploit one feature of TASEP-LK, which we state 
by saying that the binding-unbinding process can 'decouple' 
the currents at the segment boundaries. Indeed, for the com- 
posite LD-HD phase it is known [17] that the in-current J~ 
depends on the in-rate a eff only, whereas the out-current J + 
depends on the out-rate /3 eff only. Any LD-HD segment there- 
fore lifts one coupling constraint in Eqs. O, since the in/out 
currents J ± are determined locally by the junction densities 
p v and the local connectivity. 

When complete decoupling is achieved, as is expected at 
high values of Q, one can directly deduce the junction densi- 
ties for an arbitrary network. As presented in the supplement 
lElll . this leads to an exact solution of the mean field Eqs. (Q]). 
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FIG. 3. (color online) Average current as a function of the exchange 
parameter Q, for three different (average) connectivities (c = 2, 
c — 10 and c — V oo) at K — 1.5. Agreement between simu- 
lation (symbols, for L — 400) and mean-field results (solid lines, 
Eq. (1)) is excellent. The (red) dashed line is a simplified mean-field 
result (Eq. (2)) while the dotted line denotes the Langmuir current 
J/p = pi(l — pi). Results are for a single graph instance of (9(10 2 ) 
number of junctions. 



phase. Moreover, at high ft one notices a specific alignment 
of the effective rates as given by the decoupled Eqs. ©. 

In the following section we show how effective rate plots 
allow to rationalize the scale at which density heterogeneities 
appear in the network. 
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For K > 1 we have 

f C 1 - y/i^j) t f ( 2 ) 



In Fig. [3] we compare these analytical predictions based on 
complete decoupling to both a (numerical) full mean-field so- 
lution to Eqs. 0} and to simulations. For low values of Vt 
the fully coupled description Eqs. (Q} is necessary (especially 
for low mean connectivity c). Surprisingly, the decoupled 
description Eq. © is excellent even down to relatively low 
values of ft. This is an important result, as it shows that 
single segment TASEP-LK II17I1 suffices to describe transport 
through any complex network for a wide ft range. For com- 
parison we have also indicated the average mean-field current 
in the c — )> oo limit, i.e. the TASEP-LK current which is main- 
tained even when all junctions are blocked (a eff = /3 eff = 0). 
This constitutes a lower bound to any current in any network. 

Effective rate plots Here we introduce effective rate plots 
as a way to understand intuitively the physics of active trans- 
port through networks by allowing to visualize the whole sta- 
tionary transport state of the network. In Figs. 4 we map the 
effective rates (a s , f) s ) of each segment s, obtained by numer- 
ically solving of the full mean field Eqs. (1), onto the single 
segment phase diagram. Note that the scattering of effective 
rates is due to the irregularity of the networks considered. 

For sufficiently high c Figs. |4] reveal that the effective rates 
cluster close to the origin, in the LD-HD phase; this explains 
why the simplified Eqs. © work well for high c in FigJU 
When increasing ft the TASEP-LK phase diagram changes. 
In particular the LD phase reduces in favour of the LD-HD 



FIG. 4. (color online) Effective rate diagrams for irregular graphs of 
mean connectivity c{dXK — 1.5, for given values of Q) . Effective 
rates follow from solving Eqs. d). Single graph instances consisted 
of C>(10 2 ) junctions. 



Heterogeneities for TASEP-LK on irregular networks The 
parameter ft regulates the way particles distribute along the 
network, at overall density pi, and thus determines how het- 
erogeneities develop. We characterize the stationary state 
from the effective rate plots by determining the fraction of 
segments occupying the corresponding phases, see Figs. [5j 
A complementary point of view is given in Fig. [6] with the 
distribution W(p s ) of the mean segment densities p s in the 
network. From these figures we conclude that heterogeneities 
develop throughout the network in three successive regimes 
(we discuss the case K > 1 and refer to B23I1 for K < 1): 



c=2 



c=10 




FIG. 5. (color online) Fraction of segments in LD, HD and LD- 
HD phases, in the effective rate diagram of mean connectivity (a) 
c = 2 and (b) c — 10, for K — 1.5. The transition Q c between 
the network and the segment regime (denoted by the vertical dotted 
line) is determined by the condition that tild vanishes. The graph 
instances used are those of Fig. |4] 
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( i) the network regime, for low exchange rates Q, is charac- 
terized by the presence of LD and HD segments. The distribu- 
tion W(p s ) is marked by the LD and HD peaks (whereas LD- 
HD coexistence segments are distributed evenly over the in- 
termediate density range). This bimodality implies a strongly 
heterogeneous density at the network scale. 

( ii) in the segment regime, for intermediate exchange rates 
all LD segments have disappeared in favor of LD-HD seg- 
ments. The distribution W(p s ) is dominated by the LD-HD 
peak. Although all segments have similar average densities, 
the presence of domain walls implies strong inhomogeneities 
on the segment scale. 

( Hi) in the Langmuir regime, for large exchange rates f2, the 
Langmuir phase dominates. All segments have homogeneous 
densities except for small regions near the boundaries, and no 
heterogeneities arise beyond the scale of a few sites. 




Q=0 



(equilibrium for Q->oo) 



FIG. 7. Three regimes with heterogeneities on a network, segment 
and site scale, respectively, according to the exchange parameter Q 
(see main text for Q c and fi*). Particle densities are coded in grey 
scale. Pictorials on top indicate the scale at which continuity equa- 
tions are solved. 
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FIG. 6. (color online). The MF distribution W(p s ) of segment densi- 
ties for an irregular graph instance with mean connectivity (a) c = 2 
and (b) c = 10, of O(10 4 ) junctions, for K = 1.5. At low val- 
ues the bimodal distribution of 1 11] is identified. When increasing Q 
the center peak gradually grows while the edge peaks shrink. Net- 
work heterogeneities eventually disappear at intermediate Q = Q c , 
leading to a unimodal density distribution. 



The transition between the network and segment regimes 
is sharp (identified as the point ft c where all LD segments 
disappear, tt-ld = 0). Moreover this transition has an up- 
per bound = 1/2 + f(K) \n[f(K)/(l/2 + f(K))} (with 
f(K) = \K — 1\/(2(K+ 1))), for which the LD phase disap- 
pears from the one-dimensional phase diagram 112311 . In con- 
trast, the crossover from the segment regime to the Langmuir 
regime is progressive. 

Conclusions We have analyzed active transport on a dis- 
ordered network which is immersed in an infinitely diffusive 
bulk reservoir, as may be considered a simple model for cy- 
toskeletal transport by molecular motors. The dimensionless 
parameter ft characterizes the competition between active and 
passive transport, and how the associated collective effects 
lead to strong spatial inhomogeneities. 

Three regimes arise, according to the scale at which hetere- 
ogeneities appear in the network: a network, a. segment and a 
Langmuir (site-dominated) regime, see Fig. [7J Interestingly, 
these scales also set the complexity that characterizes the the- 
oretical analysis. In the network regime TASEP-LK transport 
is coupled throughout the network, whereas in the opposite 



Langmuir regime the physics is essentially determined by the 
attachment/detachment process. In the intermediate segment 
regime decoupling implies that the transport characteristics of 
each segment follow readily from those of a single segment. 

Effective rate plots allow to intuitively understand transport 
processes through networks from the single- segment trans- 
port characteristics; from the scattering of rates over both the 
LD and HD zones we can directly deduce the role of strong 
heterogeneities, see Fig. [4] This approach yields valuable a 
priori insight into yet more complex excluded volume trans- 
port such as TASEP with extended particles 112511 . TASEP with 
multiple species 12611 and, as we show in 112311 bidirectional 
motion 112711 : as in TASEP-LK, the single- segment phase dia- 
gram serves as a basis for deducing the behaviour on the net- 
work. From the effective rate diagram approach it becomes 
also clear that our results extend to types of disorder other than 
topological which are relevant to biological systems, e.g. dis- 
order in the actions of particles at the junctions 11241 128 1. An 
interesting open question is how finite diffusion 112911 could be 
handled within our approach. 

Several conclusions may be relevant for modelling cy- 
toskeletal transport. First, we have shown that strong inho- 
mogeneities in the spatial distribution of motor proteins form 
for a wide range of parameters. Even in the case of infinitely 
fast diffusion considered here they resist the equalizing effect 
of bulk diffusion. Inhomogeneities would therefore be even 
more relevant for finite diffusion. Second, the presence of 
some exchange of motors between the cytoskeleton and the 
cytoplasm may in fact simplify a theoretical description, since 
the approximation of decoupling Eqs. (2) yields excellent re- 
sults for large enough ft > Q c (see Fig. 3 and 7). Third, 
our analysis hints at a way to regulate the spatial distribution 
of motors, and therefore their cargos in the cell, by way of 
modifying the exchange parameter Q. It can be controled in 
independent ways, via uup or uoa (through the bare biochemi- 
cal rate or through the motor concentration), or via the length 
dependence in ft: regulating the cytoskeleton mesh size, for 
example by crosslinker proteins, would allow to control the 
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length scale of heterogeneities. Values reported in the litera- 
ture HQ show that the values used here (ft, K ~ (9(1)) 
is a reasonable order of magnitude, and it is therefore tempt- 
ing to speculate that a moderate regulation of ft might indeed 
allow to provoke a crossover between the various regimes in 
living cells. 
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